Two- Component Coulomb Glass in Insulators with Local Attraction 



(N 
O 

>>■ 



C/5 



I 

O 

o 



> 
O 

o 
d 



X: 



Joe Mitchell/ Anirban Gangopadhyay,^ Victor Galitski/' ^' and Markus Miiller"^' * 

^ Center for Nanophysics and Advanced Materials, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111, USA 

Joint Quantum Institute, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111, USA 
"^Kavli Institute for Theoretical Physics, University of California Santa Barbara, CA 93106-4030 
"^The Abdus Salam International Centre for Theoretical Physics, P. O. Box 586, 34151 Trieste, Italy 

Motivated by evidence of local electron-electron attraction in experiments on disordered insulat- 
ing films, we propose a new two-component Coulomb glass model that combines strong disorder and 
long-range Coulomb repulsion with the additional possibility of local pockets of a short-range inter- 
electron attraction. This model hosts a variety of interesting phenomena, in particular a crucial 
modification of the Coulomb gap previously believed to be universal. Tuning the short-range inter- 
action to be repulsive, we find non-monotonic humps in the density of states within the Coulomb 
gap. We further study variable-range hopping transport in such systems by extending the standard 
resistor network approach to include the motion of both single electrons and local pairs. In certain 
parameter regimes the competition between these two types of carriers results in a distinct peak in 
resistance as a function of the local attraction strength, which can be tuned by a magnetic field. 
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I. INTRODUCTION 

Disordered films with superconducting correlations 
host an amazing variety of interesting phenomena such as 
superconductor-insulator transitions tuned by disorder or 
an external magnetic field, with rather unusual transport 
properties^"— . These phenomena led to a lot of interest- 
ing theoretical workii"— , which were to a great extent 
spurred by an experimental feature of many strongly dis- 
ordered films with superconducting correlations - a giant 
magnetoresistance peak - for which a full theoretical un- 
derstanding is still lacking. Nevertheless, it was already 
understood in Ref. 2i that the presence of this peak, ac- 
companied by Hall measurements, suggested the survival 
of some local pairing deep in the insulator, which is only 
gradually destroyed by an increasing magnetic field. Very 
recently, the presence of localized pairs has been verified 
by more direct STM spectroscopy^- showing the absence 
of coherence peaks in the tunneling density of states de- 
spite the presence of a superconducting gap - a fact pre- 
dicted theoretically earlier-ii. The experimental evidence 
at hand support a distinct transition from a Bose insu- 
lating phase to a Fermi insulator and clearly require a 
detailed study of the strongly insulating regime which 
incorporates survival of the localized pairs. 

A prominent material exhibiting the above phe- 
nomenology is InOa;, a commercially important and ex- 
tensively studied semiconductor. Despite the uncertain- 
ties about its complex band structure, it is widely be- 
lieved that the carriers in InOx originate from oxygen 
vacancies, likely partially compensated by the triply- 
negatively-charged indium vacancies. The recent ab ini- 
tio study of Ref. [l^ calculated the formation energy of 
oxygen vacancies with different charge and found that a 
doubly-charged vacancy has the lowest formation energy 



in a crystalline environment (in zero field) . The energeti- 
cally next best state is an empty site, while a single occu- 
pied site corresponds to highly excited state. It is quite 
plausible that this tendency for local 'pair' formation un- 
derlies the superconductivity in this system, similarly to 
compounds like PbTe^, where local negative U interac- 
tions have been proposed to lead to a non-standard type 
of superconductivity of preformed hard core bosons2i. 




FIG. 1. Illustration of the two component model: The energy 
landscape is due to the combination of on-site disorder and 
Coulomb interactions. The arrows indicate typical hopping 
processes relevant for the complex low T transport in the two 
component Coulomb glass. 

A serious difficulty for the theory of such insulators 
stems from the need to treat strong disorder and electron 
pairing effects on equal footing. In addition, recent ex- 
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periments indicate that long range unscreened Coulomb 
interactions, often neglected in theoretical approaches, 
do in fact play an important role in several materials. In 
particular, the temperature dependence of the resistance 
in strong disorder^S.^ or on the high-field side of the mange- 
toresistance peak in TiN^ and InOa—^, is often well de- 
scribed by the Efros-Shklovskii law suggesting variable- 
range-hopping (VRH) in the presence of a Coulomb 
gap-- . Although local pairing attractions, often captured 
through a negative-?/ Hubbard model, have been studied 
quite extensively in previous research focussing on the 
superconductor- insulator transition 1^,24^25 ^ such studies 
have predominantly neglected long range Coulomb inter- 
actions. In contrast, we do include the latter and focus 
on more insulating regimes where Coulomb interactions 



H = X!'^*^^+2 ^ —{ni-v){nj-v)+'^ -Ui{B)ni{n, 



In the above, = 0{W) represents the random on- 
site potential due to the disorder, e.g. in the form 
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of randomly positioned dopants, f- is the unscreened 
Coulomb repulsion between the localized carriers and 
Ui{B) is the local pairing interaction renormalized by 
the Coulomb repulsion between charges localized on the 
same site (within one lattice spacing). We assume that 
it is tunable, e.g., by the magnetic field. The last two 
terms represent quantum hopping of the single electrons 
as well as of pairs of electrons. While a pair hopping 
term is obviously generated as a second order process in 
single electron hopping, the relation between the single 

electron's tunneling amplitude t'f'^ and that of the trans- 
fer of a pair may not be simple in the real materials, 
since it may involve details of the local electronic struc- 
ture, which is responsible for the negative U interaction. 
We therefore allow for a independent pair hopping am- 

f2) 

plitude ti-' . When we will discuss transport, we will take 
the two hopping amplitudes as independent phenomeno- 
logical parameters, which translate into two independent 
localization lengths for localized single electron and pair 
excitations. The magnetic field enters the hopping terms 
via the phase factors 9ij{B). 

To solve this full quantum Hamiltonian would be an 
extremely ambitious goal. We will instead isolate indi- 
vidual aspects of this complex problem. In this work, 
we make two simplifying assumptions — firstly, we fo- 
cus on the regime in the phase diagram of these films 
where the electron pairs are indeed formed locally, but 
are far from condensation. In technical terms, we treat 
the hopping terms in the Hamiltonian under the approx- 
imation ij-j'^"* ^ max(T/F, ^) (a being the lattice con- 
stant) and thus restrict ourselves to a classical model 
where transport is primarily through thermally-induced 



play a crucial role and compete in a non-trivial way with 
the local negative U attraction. Our study also has impli- 
cations on Coulomb glasses in granular materials, where 
multiple occupation of sites is allowed. These aspects 
have recently been analyzed in closely related works^^'?2,. 

The above experimental motivations lead us to intro- 
duce a lattice model that captures the various possible 
ingredients present in the actual materials: strong dis- 
order, local attraction of electrons (favored double oc- 
cupancy of sites), and long-range Coulomb interactions, 
together with quantum transport captured by intersite 
nearest-neighbor hopping. The full Hamiltonian for such 
a system can be written in a general form, with the ex- 
perimental tuning parameters (disorder strength W and 
magnetic field B) explicitly written, as follows: 



-1) + E (^S^e^'^(^'clc, e2^'^(«)clclc,c, +h.c.) (1) 

iv) 

I 

variable-range hopping, among exponentially localized 
states. This is closely analogous to the standard analysis 
of doped semiconductors^*. When discussing variable- 
range hopping transport the hopping terms are taken into 
account via the (average) localization lengths, ^1 and <^2, 
of the single electrons and pairs, respectively, which are 
a result of the i?-dependent hoppings t'^^''^\ Obviously, 
this approach prevents us from capturing superconduc- 
tivity within the model. Nevertheless, many interesting 
physical phenomena observed in experiments, such as the 
giant magnetoresistance peak, often occur rather deep in 
the insulating phase — ^ where such a strongly localized 
approach is meaningful. 

Secondly, we assume that the entire effect of magnetic 
field is to tune the local pairing interaction. It is rea- 
sonable to assume a monotonic decrease of the pairing 
strength U with increasing magnetic field. However, in 
the present work we do not include the effects of the 
B-dependence in the hopping (orbital effects) , but we fo- 
cus entirely on the effect of changing the pairing interac- 
tion on various physical observables, such as the density 
of states and longitudinal resistance. In an accompany- 
ing work^S, we study instead the magnetic field depen- 
dence introduced by the phases in the hopping terms in 
Eqn. [T] through explicit evaluation of the B-dependence 
of the localization lengths ^1.2. In reality both effects 
are present simultaneously. We find that they both con- 
tribute to a non-monotonic magnetoresistance. 



II. MODEL 

We now focus entirely on this two-component Coulomb 
glass model, which will be shown to feature a significantly 
richer variety of phenomena than the canonical Efros- 
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Shklovskii model. The latter considers a lattice of sites, 
j, with random on-site energies for electrons, (/>i, popu- 
lated with a filling factor, v. Each site i can host only 
Ui € {0, 1} electrons. The (classical) electrons repel each 
other with unscreened Coulomb interaction e^/r, and the 
disorder is assumed to be distributed over a typical range 
W , e.g., uniformly in (fn G [— W, W]. An important hall- 
mark of such systems is the soft Coulomb gap in the single 
particle density of states (DOS), p{E), close to the Fermi 
level. For many materials with compensated doping, in- 
cluding InOj;, the disorder is strong, i.e. W ^ e^/a, 
where a is the typical distance between neighboring elec- 
trons. In that case, the Coulomb gap is theoretically 
predicted^i^ and empirically foundM to be essentially 
universal at low energies: p{E) exhibits linear variation, 
p{E) = The co-efhcient a is basically independent 

of the type of lattice, the filling fraction, and the details 
of the disorde r'^^'^'^ . We find a value a « 0.35 ± 0.01 
consistent with previous numerical studies^i^, but sub- 
stantially smaller than Efros' analytical estimate 2/tt—. 
The standard Coulomb gap shows up in transport as a 
stretched exponential resistance of the form 



R{T) - Ro exp 



The co-efficient in the exponent. 



C 



(2) 



(3) 



involves just one additional parameter: the average local- 
ization length, ^1, of single particle wavefunctions, apart 
from a numerical constant, whose value 4 < C < 5 can be 
extracted from a percolation analysis of random resistor 
networks^*, as well as from Monte Carlo simulations 
Now we extend the Efros-Shklovskii model by allowing 
double occupancy and electron pairing (cf. Fig. [1] for an 
illustration) with the classical Hamiltonian 



H 



E 
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-(ni-i/)(nj-i/)+^ Yni{ni~l), 

i 

(4) 

where rii e {0,1,2}. As mentioned above, the hopping 
will be reintroduced as a small perturbation to describe 
transport later. The local attraction energies, Ui, for 
doubly occupied sites will be our control parameters driv- 
ing the crossover from the electron-dominated regime {U 
large and repulsive) to the pair-dominated regime (U 
large and attractive). In between, we find a mixture 
of gapless single electron and pair states, which exhibits 
distinctly unique features that can be captured in ex- 
periments. Note that the model is also of interest 
for semiconductors in which doubly occupied sites (the 
upper Hubbard band) play a significant role^^i^. Many 
of the effects found here generalize in modified form to 
granular systems as well. 



III. SINGLE SITE DENSITY OF STATES 
A. Definitions 

We start by analyzing the static properties of the two 
component electron glass. We consider the single site 
density of states (DOS) within typical metastable states. 
The latter are defined as classical occupancy configura- 
tions which are energetically stable with respect to moves 
of single electrons, pairs of electrons, as well as with re- 
spect to the formation of local pairs by combining two 
single electrons, or the reverse disintegration process. Let 
Sn be the set of sites with occupancy n G {0, 1, 2} in such 
a local minimum configuration. We refer to the total en- 
ergy to add (remove) a single electron on site i as 
{EI"), and as Ef^ {Ef~) for pair excitations. We define 
the DOS for electron (pair) excitations, Pm=i(2)j 

P,niE) = ^ E SiE~Er^) + ^ J2 ^^E^ET-) (5) 



where N is the number of lattice sites, Y,^ = 5*0 U 5*1, 
T,^ = 5'i U 5*2, = Sq, and = 5*2. 

In the model without double occupancy (U — >■ oo), im- 
posing stability with respect to all possible single-electron 
moves. 



E. 



1+ 



>o. 



(6) 



is sufficient to induce the Coulomb gap in the DOS. Ad- 
ditional multi-particle constraints impose weaker condi- 
tions and have been shown to not significantly affect the 
low-energy profile of the DOS. In contrast, we show be- 
low that the presence of double occupancies results in ad- 
ditional constraints, which affect the Coulomb gap very 
significantly. In the following, we describe the evolution 
of the DOS as the attraction strength is tuned. 

Clearly, for strongly repulsive U, when all double occu- 
pancies are forbidden, the system reduces to the standard 
Efros-Shklovskii model where the canonical Coulomb gap 
with slope a w 0.35 is found in the single particle density 
of states pi. 



B. Spatially Uniform Interaction - Anomalous 
Coulomb gap 

The case of a uniform pair interaction, Ui ^ U Wi, can 
be understood essentially analytically. Fig. [5] illustrates 
the corresponding evolution of the DOS's with local inter- 
action strength U, which were obtained from numerical 
simulations which we will describe below. In the attrac- 
tive case, U < 0, all electrons remain paired in local 
minima, i.e. sites are either empty or doubly occupied. 
This is so because any singly occupied site would lower 
the energy be admitting a further electron brought in 
from far away. The pair-DOS, p2{E), is linear at low E, 
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FIG. 2. DOS for difTerent uniform interaction U. (7 = is 
a critical point at which both pi, 2 have a linear pseudogap. 
The slope of the single particle DOS pi is suppressed to a/Ae^. 
For net repulsion, U > 0, pi has the canonical slope a/e'' at 
lowest energy, followed by a hump at the scale U /2, crossing 
over to the critical slope, while pairs are gapped up to i5 = U. 
For U < 0, single electrons have a hard gap |?7|/2, while pairs 
are pseudogapped with slope a/lGe**. Note: For these plots, 
the chemical potential was explicitly zeroed when averaging 
the DOS over the various initial occupancy-distributions 



with the canonical slope -jj-yc corresponding to charges 
2e. This results from the pair stability constraint analo- 
gous to Eq. dni), 



Ef+ - E]- - 4ey > 0. 



(7) 



This condition automatically ensures stability with 
respect to single-electron moves and pair forma- 
tion/disintegration and thus constitutes the dominant 
condition determining the low-energy pair-DOS. These 
assertions are easy to check case by case, using that for 
U < single particle excitations are given by 



i± 



Ef^±\U\ 



{U < 0), 



(8) 



on empty and occupied sites, respectively, being gapped 
up to energies Eg = As compared to a pair move, 
a single particle move does not only cost more in terms 
of onsite energy per particle, but also gives back less in 
terms of the polaronic interaction term Cij. Likewise, one 
checks that if (O is satisfied, it is always unfavorable to 
let a pair disintegrate into two single electrons, partly 
because one loses the attraction energy U, and partly 



because one does not gain as much polaronic energy back. 
The relation implies that the single particle DOS is 
given by 



Pi{E)^2p2{2E-sgn{E)\U\), {U<Q). 



(9) 



From this it follows that the single-DOS at energies be- 
yond the gap and close to it goes as pi{E) = -^{\E\ — 

Eg). ^ 

The point of no net interaction, U = 0, constitutes 
a critical point, where both pi{E) and P2{E) have soft 
excitations near E = Q. However, most remarkably, the 
slope of pi is reduced by a factor of 4 from its univer- 
sal value a/e^ in the canonical model, as if it were the 
Coulomb gap of a system with effective charge e* = v^e. 
This geometric mean of 2e and le arises because the gap 
imposed by the pair constraints ([7]) is probed by le ex- 
citations. Indeed, for each pair of sites admitting a pair 
move, the constraint 



E]- 



2eii > 



(10) 



must hold, as one obtains by inserting ^ ioi U — 
into ([7]). This is indeed a more stringent constraint than 
Eq. (HI). 

On the repulsive side, U > 0, pairs are gapped up to 
energy Eg = U. Mathematically, this follows simply from 
the fact that on empty sites, one has 



U, [/ > 0, 



(11) 



with E^^ > 0, and an analogous relation for doubly oc- 
cupied sites. Indeed, to accomodate a pair in a potential 
well, the well must be at least as deep as —U, which en- 
sures that the second electron is just loosely bound. The 
minimum energy required to remove the pair from such 
a well is U . Similarly, injecting a pair into an empty site 
costs at least the repulsion U of the second electron. 

On the other hand, for repulsive U, pi{E) remains 
ungapped. At low energies [\E\ <^ the universal 
single-electron Coulomb gap with slope ^ emerges: in- 
deed, the vast majority of stability constraints involving 
sites at these energies are single-electron constraints. At 
larger energies, \E\ ^ J7, one can ignore U in the sta- 
bility constraints, which then reduce again to Eqs. (I7ll0p 
and thus lead to a slope of ^ (for E below the Coulomb 

gap, i?cb ~ w >' This immediately leads to an inter- 
esting prediction: in the repulsive case, at intermediate 
energies, U/2 < \E\ < U, pi{E) is non-monotonic, as is 
indeed confirmed by the numerical data in Fig. H 

Let us now characterize the single-DOS in this regime 
in more detail. At small positive E <^U, pi{E) receives 
essentially equal contributions from empty and singly oc- 
cupied sites. Likewise, for negative energies, it receives 
equal contribution from doubly and singly occupied sites 
. If we denote the respective contributions as pf\E), 
p'^i^ {E) and pf^ (E) (with superscripts denoting occupan- 
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FIG. 3. Breakup of the single-particle DOS = + 

Pi^^(s) + Pi'^^(^) for repulsive U = 0.7, split according to the 
site occupancies, as described by Eqns. 1121 1131 and 1141 Note 
that the contribution to pi (e) from singly-occupied sites ends 
at e = ±?7, exactly where the pair-DOS /52(e) begins. Inset: 
The dip in pi(e) corresponds to p^i\() going to zero at e = 
±17, as emphasized by the dashed lines. 



cies), we find empirically that 



2 e^'- 
a \E\ 
2~ 



< < C/, (12) 



i.e., the ground state occupation is practically uncorre- 
lated with the excitation energy. 

From Eq. it follows that the pair DOS P2{E) sat- 
isfies 



(13) 



Thus, from ([T2|) . beyond the pair gap it starts off as 



P2{E) w - — -J — , 



(14) 



as can be seen in Fig. |3] 

Note that the contribution to pi [E] from singly occu- 
pied sites is restricted to the energy range \E\ < U, since 
otherwise spontaneous particle rearrangements would oc- 
cur. Further, p\^\e) satisfies the simple relation 



p\^\e) = pY>{E ~U), 0<E<U, 



.(1 



(15) 



which expresses the relationship £'^+ = E^~ + U > be- 
tween particle addition and removal for singly-occupied 
sites. Note that this implies in particular that 



p['\\E\^U) 



2e4 



(16) 



tends to zero at E — ±U, and has a maximum around 
E = ±U/2. 



At the same time the contributions pf''^\E) to the 
single-DOS do not exhibit any sharp features at energies 
of order U , except that they smoothly roll-over from a 
slope at |S| < [/ (cf. Eq. ((H])), to a slope that ap- 
proaches ^ ioi \E\ > U . As a result, the full single 
particle DOS, which is the sum of these two contribu- 
tions, exhibits a local maximum around E — ±J7/2 and 
a local minimum around E = ±t/, essentially reflecting 
the properties of p''i\e) imposed by the extra constraint 
psp . Very similar physics was found recently in granular 
systems — where the occupancy of sites is nearly un- 
limited, in which case Eq. ([15]) applies essentially to the 
whole 1-particle DOS, and imposes mirrored Coulomb 
gaps. 

Despite the absence of quantum fluctuations, the de- 
scribed evolution of the DOS has a lot in common with 
quantum critical phenomena — , where U plays the role 
of the detuning parameter from criticality. The critical 
behavior is restored at energies \E\ ^ |C/|, with linear 
DOS-s and anomalous slope of pi{E). At low energies, 
the non-critical phase appears, where one type of carriers 
is gapped out, while the other type exhibits a universal 
Coulomb gap. We also note that the features of the DOS 
and the underlying mechanisms found here have similar- 
ities with those in a recently proposed model of strongly 
and weakly interacting two-level systems^. 



C. Numerical simulations 

In order to analyze further details of the DOS, as 
well as the case of random local interaction Ui, we per- 
formed numerical studies. To study metastable states, 
we start from a random configuration of occupancies, rii 
(€ {0; 1, 2}), on a half-filled triangular lattice of size 200 
X 200. We choose a triangular lattice with commensurate 
filling so as not to introduce extra strain in the system in 
the limit of weak disorder (note that if the filling is not 
commensurate there is still some strain from the lattice) . 
However, we focus on strong disorder where the effect of 
the lattice type is expected to be small. 

Following a similar protocol as described in^, we al- 
low re-distribution of occupancies through single particle 
moves, pair moves and pair dissociation/formation — the 
last one within a restricted spatial range — that lower the 
total energy of the system until the system stabilizes in a 
local minimum of the energy. In this context, it is impor- 
tant to recall that the appearance of the Coulomb gap in 
the single particle DOS does not require stability with 
respect to multi-particle moves, and is not very sensitive 
to the latter. This is because the single-particle moves 
impose the strongest stability constraints'^'*. By a similar 
reasoning, the universal features in the DOS for uniform 
U in the model considered here result from single particle 
and pair stability constraints. It is thus reasonable to ex- 
pect that the class of moves considered above imposes the 
strongest stability conditions, determining the essential 
features of the single site DOS-s, pi{E) and P2{E), for 
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single electron and pair excitations, respectively. Further 
multiparticle processes may relax the system to lower ly- 
ing metastable states; however, such states are expected 
to have very similar single site density of states and trans- 
port properties. The single site DOS was obtained by 
calculating the histogram of the energies to add or sub- 
tract an electron or pair from each site, cf. Eq. ([S]). These 
DOS-s were averaged over many different disorder real- 
izations, typically of the order of 100 for the 200 x 200 
sized systems. 

We measure all distances in units of the lattice constant 
a, and in our finite-sized samples, the intersite distance 
r,j , has been chosen as the minimum distance on a torus 
defined by periodic boundary conditions. Energies are 
measured with reference to the chemical potential fi in 

2 

units of the nearest neighbor Coulomb repulsion ^ . The 
chemical potential in this case is determined as the av- 
erage of the smallest energy to add and remove an extra 
particle from a given metastable state. Pair energies are 
measured from the reference energy 2/i. 

We choose the on-site disorder to be randomly dis- 
tributed in the interval [— W, VF]. It is well-known that 
a disorder of order unity or more is required for the 
DOS to tend to an essentially universal Coulomb gap, 
pi{E) = at low energies. In our model, since a site 
is allowed to have double occupancy, the strong disor- 
der condition is met when W exceeds the typical nearest 
neighbor interaction of two doubly occupied sites. In 
order to find DOS features which approach a universal 
limit, we therefore chose to work with disorder W = A. 
At substantially weaker disorder the low energy DOS's 
were indeed found to be non-universal. 



D. Spatially Disordered Interaction 

In a disordered system, it is more realistic that the 
pairing energies Ui are non-uniform. Fig. H] shows the 
result for the DOS-s for a model with Ui distributed 
randomly in an energy range [U - l^U,U + l^U]. The 
sharp features of Fig. [5] are smoothened. The gaps in 
Pi{E) and p2{E), for repulsive and attractive U respec- 
tively, are smeared out. Low energy single particle states 
leak into the gap of pi{E) as soon as there are positive 
C/j-s. The density of such states grows with increasing 
U, and eventually saturates to the standard linear pseu- 
dogap with slope Closely analogous considerations 
apply to P2{E) upon decreasing U. The detailed behav- 
ior in the intermediate regime \U\ < AU is complicated 
and presumably non-universal. 

The non-monotonic humps in pi discussed in sec- 
tion 11116] survive only if AU is sufficiently small as com- 
pared to U. Thus they are probably best sought after 
in crystalline samples, where the local environment of 
different impurities are similar, giving rise to a narrow 
scatter, AU. 




Energy Energy 

FIG. 4. Comparison of the DOS's for constant (left) and 
random (right) U, with strong scatter AU — 2 (units of e^/a). 
The sharp gaps and humps are smoothed out by disorder, but 
the overall trend of increase/decrease of the low energy DOS 
remain intact. 



IV. TRANSPORT AND RESISTANCE 

A. Choice of parameters 

In an insulator, a reduced density of states is usually re- 
flected in an exponentially increased resistance. It is thus 
interesting to ask what happens in the "mixed regime" 
of our model, where both pair and single electron excita- 
tions are ungapped. If transport was dominated by one 
type of excitation only, one would expect an increase of 
resistance upon approaching the mixed regime from ei- 
ther side, since the DOS of the dominant carrier type 
diminishes. However, transport is more complicated in 
this two-component Coulomb glass. Electron and pair 
hops do not take place independently in the sample, but 
combine to form a network of interconnected pair and 
single electron moves, as illustrated in Fig. [TJ Trans- 
port is a complex functional of the combined density of 
states. In order to study this insulating regime, we have 
generalized the construction of an effective network of 
Miller- Abrahams resistors^ which include both pair and 
single particle processes. In order to elucidate the inter- 
play between pair and single particle transport, we have 
neglected spin blocking effects in the random resistor 
network'^^ (as may be justified in strong spin orbit cou- 
pled materials). Spin dependent effects may be consid- 
ered elsewhere. The elementary hopping resistances were 
evaluated in a mean field fashion for selected metastable 
states2^. 

Details of the resistor network, the required steps and 
approximations, are described in the appendix. The only 
important point to note is the exponential dependence of 
the effective resistances of the network on temperature 
and localization length. This restricts the accessed range 
of energies and typical hopping distances of the electrons 
participating in the network and using a percolation ar- 
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gument^^, one can determine the functional dependence 
of the resistance on temperature. 

In presence of the Coulomb gap in the density of states, 
and if one type of carriers dominates the low T transport, 
the functional dependence is of Efros-Shklovskii type, cf. 
Eq. ([2]), with 



Tn*'^ = C 



(17) 



where (5i=i.2 is the charge of the carriers in units of e, S.i 
their average localization length and C « 4 — 5. These 
localization lengths may in principle be evaluated from 
an analysis of the elementary localized excitations above 
the ground state, whose spatial extent is governed by 
the hopping terms in Eq. [T] The magnetic field enters 
this localization length via phases in the hopping and re- 
sulting interference effects, as discussed e.g. in Refi^i^. 
However, for the purpose of analyzing the effect of vary- 
ing pairing strength, we assume the localization lengths 
to be constant. This may describe very well an experi- 
mental situation in which the local interaction U is tuned 
(by chemical modifications or gating), without affecting 
the localization lengths. In contrast, in the case where U 
is tuned by a magnetic field, the effects described below 
will necessarily be superposed over quantum interference 
effects, which affect rather sensitively and may well 
dominate the effects which we address below. 

As we shall explain below, under certain circumstances 
we obtain a nonmonotonicity in resistance as a function 
of U. The latter is most prominent when we have |^ = 4, 

(1 2) 

in which case the Efros-Shklovskii temperatures Tq ' 
are the same for both single-electron and pair-transport 
resulting in an interesting competition when U is tuned 
across zero. In reality the ratio varies greatly across 
the phase diagram of disordered films with supercon- 
ducting correlations. Indeed the localization length of 
preformed pairs must diverge at the transition to a su- 
perconductor, while ^1 remains non-critical^^^. On the 
other hand, far in the insulating regime, ^2 is expected 
to become shorter than ^1, because pair tunneling is sup- 
pressed. Therefore a regime where ^2 > should cer- 
tainly exist, and below we consider the particularly in- 
teresting case ^2/^1 = 4. We should keep in mind how- 
ever, that such a large ratio presumably implies relatively 
strong quantum fluctuations, due to rather important 
hopping terms. With this caveat in mind, the essentially 
classical description of the two-component Coulomb glass 
presented here should be taken as a phenomenological ap- 
proach to capture Coulomb frustration effects on a sys- 
tem with variable range hopping transport of competing 
carriers. 



B. Results 

In the case of uniform U, there is no genuine mixed 
transport regime at low T, since one of the two car- 
rier types is always gapped. Here we find the resistance, 



R{U;T), to be a flat, i.e., constant function of U within 
error bars. Despite a suppression of pi{E) (see Fig. [2]), 
the resistance does not increase significantly. As U is 
decreased, the pairs "fill in" the resulting gap left in 
the transport channels, keeping the resistance essentially 
constant. 

However, the situation is more interesting with ran- 
dom J7i, where a genuine mixed two-component carrier 
regime exists. If pairs and single-electron excitations 
have strongly disparate localization lengths, R{U) is still 
monotonic under tuning of U , essentially refiecting the 
evolution of the DOS of the less localized carrier type. 
However, the two carriers do compete significantly in an 
intermediate regime. Indeed, we find an interesting non- 
monotonicity in the resistance; see the top panel in Fig. 
[5l At low temperatures we find a relatively significant 
peak in the resistance, centered around a small U < Q. 
This feature is even slightly enhanced by increasing the 
randomness AC/. Transport in the peak region is par- 
tially by pairs, which break up and propagate as single 
electrons, and then recombine again. The maximum of 
resistance occurs when roughly an equal number of single 
and pair hops form the critical links of the percolation 
network, see Fig. O 

A plausible qualitative explanation for the numerically 
observed peak is the following: it is difficult to connect re- 
gions in which pair or single electron transport is favored, 
as opposed to the regimes \U\ >• 1, where transport is 
dominated by one type of carrier only. In other words, 
the mixed regime suffers from "contact resistances" be- 
tween pair- and single-dominated parts of the resistor 
network in the following sense. A piece of transport path 
of pairs, must be connected to two good strands of single 
particle transport, and likewise a single particle transport 
path must find another one to continue as a pair path. 
Both links require some matching which tends to increase 
the overall resistance. This phenomenological explana- 
tion of our numerical observations bears some resem- 
blance with the idea that superconducting islands may 
act as weak links in a single-electron-dominated trans- 
port regime on the insulating side of the SIT, as proposed 

™25 



V. CONCLUSION 

The model discussed in this paper, should be realized 
in disordered materials with a strong tendency for lo- 
cal attractions (negative Hubbard U). It also predicts 
interesting effects for cases where local interactions are 
moderately repulsive, such that multiple occupancy of 
sites is still possible. The occurrence of negative U inter- 
actions is likely to be concomitant with a bosonic type of 
superconductor-to-insulator transition upon further re- 
duction of the disorder. In such samples, U may also be 
tuned by an external magnetic field which has a depair- 
ing effect on the electrons. However, since a magnetic 
field also sensitively affects localization lengths, it would 
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FIG. 5. Top panel: Peak in the resistance upon tuning U 
a,t T — 0.04 vs. the average interaction U (with ^2 = 4^1 ); 
and the fraction of pair hops in the percolating cluster (with 
resistances within 30% of the percolating resistance). The 
resistance peaks when roughly half of the critical resistors are 
pair and single moves, resp. Lower panel: fraction of paired 
electrons in a typical metastable state. This fraction smoothly 
decreases across the "mixed regime", since the bulk of such 
pairs is inactive in transport. 



be desirable to use other, non-magnetic means to influ- 
ence the local interactions, too (such as pressure, chem- 
ical doping etc). If the disorder in the local U-s is large 
we find a regime around C/ ~ 0, where both pairs and sin- 
gle electrons contribute to the activated transport, and a 
non-monotonic resistance as a function of U results. 

For the strongly localized, classical limit of the two- 
component Coulomb glasses we found several interesting 
effects on the low energy density of states. In particular, 
we find that the tendency for local attraction leads to a 
suppression of the density beyond the standard Coulomb 
gap. At the point where local attraction and repulsion 
balance to produce vanishing net interaction U = we 
find that the 2d Coulomb gap is reduced by a factor of 
4 from its canonical value. More generally, if multiple 
charging of the same site with M charges (without paying 
additional local charging energy) were allowed, one would 
find a suppression by a factor of M^. 

For the case of moderately repulsive interactions U > 
0, if the randomness in the interaction energy AU is small 
compared to the average interaction U, our model pre- 
dicts the existence of non-monotonic humps in the sin- 
gle particle density of states. If pair transport is sup- 
pressed due to strong localization, the non-monotonicity 
of the single-electron DOS in the repulsive case U > 
should show up as a kink in the resistance R{T) around a 
temperature T, (?7/2)^/(Ce^/^), where it crosses over 
from an Efros-Shklovskii law with a higher value of To 
to a less steep R(T) and a twice smaller Tq at lower T. 
The humps in pi should also leave traces in AC mea- 
surements^*, or more direct measurements of the DOS 
such as photoemission or tunneling from a broad junc- 
tion^. These DOS features may also be relevant for the 



more involved experiments of memory effects in deep in- 
sulators^Si^, where doubly occupied sites with repulsive 
interactions are known to be presenlj^^. 

A measurement of the pair-DOS, especially on the at- 
tractive side U < could be attempted through mea- 
surement of the tunnelling conductance from a (wide) 
superconducting probe, similar to the experiments per- 
formed by Dynes et al^^. 

In this work we have taken the localization lengths 
to be independent of the tuning of the local interaction 
strength. If the latter is tuned by magnetic field, a full de- 
scription needs to take such quantum effects into account, 
however. In fact, it has been argued^i^ that the field 
dependence of localization lengths of pairs and electrons 
are opposite, which is probably an important ingredient 
for a strong magnetoresistance peak. Here we show that, 
on top of that effect, the complex energetics and trans- 
port phenomena in the two component Coulomb glass 
can even enhance such a peak. A discussion of the com- 
bination of both magnetic field effects (tuning U and the 
localization) is left for future studies. 
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APPENDIX A : RESISTOR NETWORK 
MAPPING 

The problem of finding the conductivity of a sample, 
where transport takes place via variable range hopping of 
electrons from singly occupied sites to empty sites, can be 
mapped to a n eq uivalent random resistor network prob- 
lem (see Ref. |4l| for a derivation) . The sites in the hop- 
ping problem can be mapped to the vertices of this net- 
work, and the inter-site transition rates to the resistances 
linking these vertices. Finding the sample resistance sim- 
ply reduces to calculating the effective resistance of the 
resistor network through a percolation approach— . 

It is possible to formulate a similar resistor network 
mapping for variable range hopping on a lattice where 
double occupancy is also allowed. Early efforts in this 
direction were made by Kamimura et. alj^. There, the 
inter-site transition rate was taken as the sum of tran- 
sition rates of the four types of possible single-particle 
processes characterized by the occupancies of the initial 
and final sites before the hop. 

In our two-component model, we additionally intro- 
duce the pair-hopping channel between two sites along- 
side the single-particle transport channels considered in 
Ref. Is^: in what follows, we refer to these single par- 
ticle and pair hops as first order processes. However, 
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with this introduction, two-particle processes involving 
more than two sites may also become significant as they 
can potentially provide lower-resistance alternatives to 
two-site pair hops between given initial and final occu- 
pancy configurations. For brevity, we refer to these two- 
particle processes involving more than two-sites as second 
order processes with the idea that two such processes in 
succession (a pair breaking followed by re-joining of the 
lone electrons) can constitute a first order process (a pair 
hop). The book-keeping of the different allowed hopping 
processes based on occupancies prescribed in Ref. |33 be- 
comes increasingly more cumbersome, as one takes into 
account these second-order processes. 

Hence we formulate a more general prescription for 
the resistor-network mapping which can be naturally ex- 
tended to include higher-order processes. We first de- 
scribe this reformulation for the case where transport 
takes place through first-order processes only. In this 
case, each vertex of the equivalent resistor network corre- 
sponds to a node {i,ni) defined as a site i together with its 
occupancy > 1 . In a system with double occupancies 
allowed each site gives rise to two nodes, while upon elim- 
inating double occupancies the only nodes are («, = 1), 
which reduce to the standard Miller- Abrahams network. 
The various first order processes between the two given 
sites i and j correspond to resistances between different 
nodes of the network: see Table |T] for a full description of 
the four possibilites. In networks without double occupa- 
tion, there is only one resistance, 1; j, 1) associated 
to a given pair of sites. 

In defining the resistances, it must be noted that, while 
the occupancies of the two sites under consideration as- 
sume all possible nonzero values, the occupancies of all 
other sites are frozen at a pseudoground state. 

The expression for the transition rate given below in 
Eq. (jl8p for a single-particle hopping process of one type 
tallies exactly with that of Ref. 37. Since it is expected 
that, between two given sites, one type of single-hopping 
process is energetically favoured, and thus has least re- 
sistance compared to the other three, the two resistor 
network mappings might be naively expected to produce 
identical results if only single-particle hops are consid- 
ered. However, a conducting path through the sample 
constructed from these least resistance hops might in- 
volve an occupancy mismatch between resistors sharing 
the same site: the alternative paradigm for the resistor- 
network construction in our approach, which includes oc- 
cupancy in the definition of vertices of the network, ex- 
cludes this possiblity. 

In the second section, we will include second-order pro- 
cesses by a simple extension of the definition of nodes 
from (i^rii) to {i,ni,j^nj). Herein lies the utility of our 
approach : higher-order multi-particle processes can be 
easily incorporated simply by extending the dimensions 
of the node-network. The resistor network constructed 
from these generalized nodes contain the first-order pro- 
cesses as a subnetwork as will be described in detail later. 
This extension, however, makes calculation of the effec- 



tive resistance very expensive. Therefore, in the last sec- 
tion, we describe possible approximations so that we can 
investigate the effect of pair breaking/formation (crucial 
two-particle moves in the regime where both single and 
pair DOS's are ungapped) despite staying within a sim- 
pler resistor- network with nodes of the form (i, n^). 



1. Resistor-network construction for first order 
processes 

Let us now describe the resistor network mapping in 
detail. We start by reaching the pseudo-ground state, 
used earlier to extract the density of states. 

In zero field, the time-averaged rate of transfer of elec- 
trons through single- hops from node («, rii) to node (j, rij) 
(note: is the occupancy of site i before the hop, while 
Uj is the occupancy of site j after the hop) is given by 

T {i,n^;j,nj) = e « P{i,ni; j,nj) (18) 

Here rij = \rj — r^j is the distance between the sites 
and ^ is the localization length of the electronic wave- 
functions for single-hops and that of the pairs for pair 
hops. E(ni , rij — 1) is the total energy of the system when 
sites i and j have occupation numbers rii and rij — 1, 
respectively (for all sites other than i and j, we use 
the pseudo-ground state occupancies). Thus, the term 
e~^^ /Z acts as a Boltzmann probability for the initial 
occupancy-configuration of the two sites. P{i,ni; j,nj) 
is the amplitude for phonon emission or absorption, as 
the electron system changes from a configuration with 
(j,ni), (j, rij — 1) to one with {i,ni — 1), {j,nj). Up 
to pre-exponential factors, which we approximate by a 
uniform value here (set to 1 by a choice of unit of time), 
P{i,nfJ,nj) is givenhy N{i,n^,j,nj) = ^g|Ag(i.4;j.^j)i _i 
(absorption, AE > 0) and (1 + N{i,rn, j.rij)) for emis- 
sion, AE < 0. The two cases can be combined into a 
single expression for N{i, rii] j, Uj): 



P{i,n,;j, n.j) 



1 



gl3AE{i,ni;j, rij) _ I 



(A2) 



Finally, Z is the "two-site partition function" ex- 
pressed as 



E 



where the sum over rii and Uj runs from to 2. 

A pair-hop from site i to site j corresponds to a resistor 
between nodes {i,ni — 2) and (j, — 2) defined by 
analogy to Eqn. [18] as 



r°(i,ni;j,nj) = e" 



e 



z 



-P{i,nf,j,n.j) 



(A3) 
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Initial occupancies 


Resistance R{i,ni; j,nj) 


Site i 


Site j 


1 





R{i,l 




1 


1 


R{i,l 


i,2) 


2 





R{i,2 




2 


1 


R{i,2 


J, 2) 



TABLE I. Description of possible single-particle hopping pro- 
cesses between two sites i and j through resistances con- 
structed from nodes of the form {i,ni) : it is important to 
note that rii is the occupancy of site i before the hop, while 
rij is the occupancy of site j after the hop 



where the variables are the same as defined in Eqn. 1181 
except that E{ni ,nj — 2) is the total energy of the system 
when sites i and j have occupancies Ui — 2 and n,j — 2 — 
0, respectively. The nodes connected by this pair hop, 
namely — 2) and (j, = 2), are also connected 

by a single hop (the last hop described in Table U and 
we choose the smaller of the two resistance values as the 
effective resistance between the vertices. 

Note the difference of Eqn. [TS] from a similar expres- 
sion for the time-averaged rate of transfer given in Ref. 
m. In the latter, for a single hop from j to j, the probabil- 
ity of occupancy = 1, Uj = is implemented through 
a product of the individual occupancy probabilities as 
/i(l — fj), where fi is the Fermi distribution function, 
instead of the Boltzmann probability-term e~^^/Z in- 
cluded here. In such an approach, detailed balance in 
the absence of an electric field can be obtained only by 
dropping the "polaron term" ^ (particle-hole interac- 
tion) when calculating AE{i,ni; j,nj) . This simplifica- 
tion leads to a somewhat different value of the Efros- 
Shklovskii temperature Tq (cf. Eqn. [TT]) as compared to 
the treatment used here^. However, the particle-hole in- 
teraction may play a more vital role in our model, where 
we include electron pairs: the increased importance arises 

2 

from the fact that the polaron term is equal to — for a 

rij 

single hop but 4 times that for a pair hop and thus may 
be rather significant in the mixed regime favoring pair 
transport as a whole over pair breaking. 

As a check of our prescription, it can be easily ver- 
ified that upon barring the polaron term, in the limit 



Here, rtj-ki is taken as min [r^^ + rji,ru + rjk], keeping in 
mind the indistinguishability of the hops {i,ni,j^nj) — 
{k,nk,l,ni) and {i,n,,j,nj) -> {l,ni,k,nk) (see figureE]). 

Although such a two-particle move potentially in- 
cludes formation or breaking of pairs, we choose for ^ 
the localization length of the single-electron wave func- 



of large U where pair formation/transport is hindered, 
the expression for the transition rate given in Eqn. [TH] 
completely agrees with the resistor network construction 
used in Ref. M (see also Ref. 

In the presence of a weak electric field, one can asso- 
ciate a resistor between the nodes (i, ri;) and {j, rij) with 
resistance value given by 

kT 

R{i,n,;j,nj) = R{j,nj;i,n,) = ■. -, (A4) 

which is guaranteed to be nondirectional due to detailed 
balance. 

2. Extension to include second order processes 

The resistor network mapping stated above allows a 
generalization to include second order processes. Below, 
we shall describe an extension to include generic two- 
particle hops and then make approximations to a specific 
set of two particle hops: namely, ones that involve the 
formation or distintegration of a pair. 




j 1 j 1 



FIG. 6. One activated resistance, R{i,rn, j,nj;k,nk,l,ni), 
corresponds to the above two hops, as the particles are indis- 
tinguishable from each other 

To describe a second-order process, we extend the no- 
tion of a node: now, each node is defined by a pair of sites 
with corresponding occupancies as (i, rti, j, Uj), subject to 
the constraint rii + nj > 2. For a two-particle move from 
node {i,ni, j,nj) to {k,nk,l,ni) (note: rii and rtj are the 
occupancies of sites i and j before the two-particle hop 
has taken place, while rik and n; are the occupancies of 
sites k and I after), the time-averaged rate of transfer of 
electrons can be given by the equation 



(A5) 

I 

tion in all cases except for pair moves, {i, rii, i, rii) 
{k,nk,k,nk) (see Ref. |H for a discussion of such is- 
sues in the context of the Hubbard model). Other 
than pair hops, the remaining first order processes are 
treated as a subnetwork of this network of generalized 
nodes in the following sense — a single particle-hop 



T"{i,ni,j,nj;k,nk,l,ni) = e 



-l)E{ni,nj ,nk — l,ni — l) 



-P{i,niJ,nj;k,nk,l,ni) 
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from site {i,ni) to {j,nj) involves a "trace" over re- 
sistances of the form R{i,ni,k,nk', j,n,j,k,nf^) with the 
indices k and running over all sites and occupan- 
cies respectively. In practice, in the percolation ap- 
proach, this "trace" is performed by "activating" resis- 
tances R{i, rii, k, nk]j, nj, fc, Uk) for all k and nk simulta- 
neously for the single-particle hop {i,ni) — > {j,nj) since 
these resistances all have the same magnitude and are 
distinguished only by their location in node-space. 

The probability for absorbing (or emitting) a phonon 
with the required energy is in analogy to Eq. (jA2p : 



P{i,ni,j,nj;k,nk,l,ni) 



1 



and the resistance associated with the link between the 
nodes is : 



R{i,ni,j,nj;k,nk,l,ni) 



kT 



e'^r'^{i,ni,j,nj;k,nk,l,ni) 



(A6) 



3. Simplified algorithms to include pair 
breaking/ formation 

The percolation algorithm with generalized nodes 
(i,ni,j,nj) takes into account all possible second or- 
der processes whose number grows roughly as the fourth 
power of the system size. This constrains the approach to 
small system sizes (upto a 50 X 50 lattice). Nevertheless, 
performing the full percolation analysis with the general- 
ized node-network for such small sizes allows comparison 
with certain approximate networks that we describe be- 
low. Once the results are seen to agree well, these latter 
approximate networks can then be used on larger systems 
to calculate the effective network-resistance through the 
percolation approach^. 

The first simplification comes from the sparsity of 
the generalized resistor network in the sense that low- 
resistance second order processes (within percolation 
threshold) involve only small-range hops. As a result, 
for instance, for a 50 X 50 lattice, while calculating the 
resistances for second order processes, we can restrict the 
range of hopping to within 5 sites only. 

Since, for reasons discussed above, we intend to 
focus on pair formation/disintegration only, another 



approximation is to retain only those resistances 
that correspond to this class of second-order pro- 
cesses. In the node-language, these resistances are of 
the form R{i,ni,j,nj;k,nk,k,nk) (pair- formation) and 
R(k,nk, k,nk',i,ni, j,nj) (pair-breaking). This imposes 
storage requirements which go as a third power of the 
system size. 

Moreover, with this simplification, we can revert to us- 
ing nodes of the form {i,ni) for first order processes. In 
effect, we are thus using a mixed definition of nodes, de- 
pending on the order of the process, which is unsuitable 
for a textbook percolation analysis. Let us therefore de- 
scribe the percolation approach in slightly more technical 
detail. We activate the resistances in increasing order of 
magnitude, irrespective of the order of the process (and 
thus the definition of the nodes connected). 

During the percolation analysis ("activating" resis- 
tances in ascending order of magnitude till a percolating 
cluster is obtained), we use the stored second-order resis- 
tances to look for an effective 'short cut' pair transport 
in the following way. While activating, say a resistance 
corresponding to pair breaking originating from site i, 
we check if the lone electrons from this pair breaking 
are connected through a path of already activated first 
order single-particle processes to another already acti- 
vated second-order pair-formation at some other site j. 
If such a path exists, we refer to it as a 'short cut' pair- 
transport. In this case, we treat the nodes (?,rti) and 
(j, Tij) as if connected through a pair-hop with the equiv- 
alent resistance equal to the last-resistance activated in 
this 'short-cut' path. The advantage of this approximate 
approach is that the percolation criterion is being applied 
effectively only to the simpler node-network with nodes 
of the form {i,ni) and thus involves a drastic reduction 
in computation time. 

With these two simplifications, it is feasible to study 
transport for system sizes as high as 200 X 200. This 
approximate algorithm has been used to check the ro- 
bustness of the results obtained in the paper to inclusion 
of the higher-order processes. We find that as antici- 
pated, the resistance for intermediate U (between single 
and pair-dominated regimes) does decrease due to the 
extra channels thus included. However, the results de- 
scribed in the main text — notably the nonmonotonicity 
of the resistance — remain present. 
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